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Spatially homogeneous cosmological models with a posi- 
tive cosmological constant are investigated, using dynamical 
systems methods. We focus on the future evolution of these 
models. In particular, we address the question whether there 
are models within this class that are de Sitter-like in the fu- 
ture, but are tilted. 



I. INTRODUCTION 

This paper is concerned with the evolution of spa- 
tially homogeneous cosmological models with a cosmolog- 
ical constant. Several authors have addressed this prob- 
lem. Stabell & Refsdal [0 investigated the Friedmann- 
Lemaitre dust models, and a generalization of their re- 
sults to general equation of state appear in [|||] . The ex- 
amination of Bianchi models with two-component fluids 
by Coley & Wainwright [Q contain the cosmological con- 
stant as a special case. The future evolution of Bianchi 
models has been considered by Wald |J , who showed that 
all non-type-IX Bianchi models with a positive cosmolog- 
ical constant isotropize. The result applies to tilted mod- 
els, but it should be pointed out that the isotropization 
is with respect to the congruence normal to the homo- 
geneous symmetry surfaces - not the fluid congruence. 
The isotropization of inhomogeneous models has been 
studied in a Newtonian context in ||, where some exact 
results were obtained. The key issue here is the degree 
to which such a constant (or effective constant, due to 
a scalar field) can lead to isotropization of the universe, 
and hence explain the presently observed near-isotropy. 
Isotropization can happen to some degree in anisotropic 
models without a cosmological constant (?]||, but not 
sufficiently to convincingly explain the degree of isotropy 
observed. Wald's theorem leads credence to the claim 
that inflationary models can isotropize the universe as 
desired, but there are limits to what inflation can achieve 
||. Hence even given all these important contributions, 
we think it is useful to examine the issue further, specif- 
ically by giving a treatment of these models in line with 
the dynamical systems methods presented in e.g. the 
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book edited by Wainwright & Ellis [juj . We will only ex- 
amine the effect of a genuine cosmological constant here; 
thus we do not consider proposals for a 'varying cosmo- 
logical constant' or scalar field (or other essentially equiv- 
alent proposals). However, we believe that examination 
of those cases by the same dynamical systems methods, 
already initiated in the case of Friedmann-Lcmaitre ge- 
ometries and more general models [pT|-^6[, will 
be a worthwhile extension of what is presented here. 

With the conventions we use, the field equations take 
the form G a b + A g a b = T a b, where A is the cosmolog- 
ical constant. In this paper, it is restricted to be non- 
negative, A > 0, with the focus on understanding the 
evolution when A > 0. Interpreting A as a vacuum en- 
ergy, this implies that the vacuum energy is never nega- 
tive. We will additionally assume a perfect-fluid matter 
source with p — (7 — l)/i as equation of state, where /i 
is the energy density, p is the pressure, and 7 is a con- 
stant. Causality then requires 7 to be in the interval 
< 7 < 2. Furthermore, 7 = just corresponds to a 
second cosmological constant and will not be considered. 

The outline of the paper is the following. In Sec. ||, 
the Friedmann-Lemaitre models are investigated. The 
detailed treatment of these well-known models is moti- 
vated by the fact that they contain many of the features 
exhibited by the more general models to be considered, 
and also illustrates the methods used in an example fa- 
miliar to most readers. Sec. [II is concerned with the 
Bianchi models. After some general statements about 
these models, we examine the simplest cases - orthogo- 
nal models of type I and II - in more detail. The section 
is concluded with an investigation of tilted LRS type V 
models, which turns out to give us reason to be cautious 
about the implications of the important Wald theorem 
In Sec. |y]we consider the Kantowski-Sachs models, 
which turn out to have a rather rich solution structure, 
with a particular interesting anisotropic boundary solu- 
tion. We end with some conclusions in Sec. [y|. 



II. FRIEDMANN-LEMAiTRE MODELS 

In order to introduce the notation used and com- 
pare with earlier results, we start by investigating the 
Friedmann-Lemaitre models. These models are homoge- 
neous and isotropic about every point, corresponding to 
a six-dimensional isometry group Gq. For a perfect fluid, 
the fluid motion then necessarily coincides with the con- 
gruence normal to the symmetry surfaces, and the models 
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are fully specified by the length scale factor S(t) and a 
parameter k G { — 1,0, 1}. The line element can be writ- 
ten 



ds 2 = -dt 2 + S 2 (t)dl k 2 , 



(2.1) 



where dl 2 is a constant-curvature three-geometry. In 
what follows, a dot denotes the derivative with respect 
to t. Defining the Hubble scalar by H = S/S, the field 
equations and equations of motion become |17|| 



The Friedmann equation: 

H 2 = i p - l 3 R + I A 
3' 6 3 

The Raychaudhuri equation: 
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The energy conservation equation: 



(2.2) 



(2.3) 



(2.4) 



Here, 3 R = 6k / S 2 is the three-curvature of the symmetry 
surfaces and q is the deceleration parameter. Assuming 
H =/= 0, we proceed by defining dimcnsionlcss variables 
according toFI 
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(2.5) 



6iJ 2 ' "" 3H 2 ' " MJY 3H 2 
The density parameter is related to K and Q\ by 
The Friedmann equation: 

n = i + K-n A , (2.6) 

while the deceleration parameter is given by 
The Raychaudhuri equation: 
37- 2 
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(2.7) 



The weak energy condition together with A > immedi- 
ately give 



< O, -1 < K, < fi A - 



(2.8) 



In addition, for models with non-positive spatial curva- 
ture 3 R, these quantities are compact, i.e. the range is 
contained in a compact interval: 



1 In [juj, the curvature variable K was defined with opposite 
sign. 



0<fi<l, -1<K<0, 0<O A <l- 



(2.9) 



Upon introducing a new dimensionless time variable de- 
fined by ' = H^ 1 d/dt, we obtain 



H' = -(l + q)H, 

K' = 2qK, 
Qa' = 2(l + g)fi A - 



(2.10) 

(2.11) 
(2.12) 



Note that the H 1 equation is decoupled, so that we obtain 
a two-dimensional reduced dynamical system for K and 
In addition, we obtain an auxiliary equation from 



The energy conservation equation: 

Q' = [2q - (37 - 2)] n. 



(2.13) 



which then shows that the Friedmann equation is an in- 
tegral of the Rayc hau dhuri equation (expl icitly , the time 
derivati ve of Eq. (2.6) is the same as Eq. ( |2.13 ) in virtue 
ofEq. 



A. Analysis of the dynamical system 



The reduced dynamical system (2.11- 2~12| ) has a num- 
ber of invariant submanifolds: 

Q = the vacuum boundary 
K = the flat submanifold 
$1a = the A — submanifold 

Despite the appearances, H — is not an invariant sub- 
manifold (As the definitio ns of the dimcnsionlcss vari- 
ables assume H ^ 0, Eqs. ( 2.10] - 2.13 ) are not valid when 
H = 0). Equilibrium points with finite values of K and 
Q\ lie at the intersection of these invariant submanifolds; 
they are 



K Q, f 



F the flat Friedmann solution 
M the Milne solution 
dS the de Sitter solution 
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Analyzing the stability of these equilibrium points, we 
find 

Eigenvalues Stability 

0<7<2/3 2/3<7<2 



F 

M 
dS 
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saddle 
source 
sink 



source 
saddle 
sink 



In addition, there might be equilibrium points associated 
with infinite values of K and/or fl A . This illustrates 
one advantage of dynamical systems with compact state 
space: there is no 'infinity' that can be difficult to ana- 
lyze. In the present case, the part of state space contain- 
ing models with non-positive spatial curvature (K < 0) 
is compact, but the region with K > is not. 
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In order to compactify this region as well, we note that 
the Friedmann equation ( |2.2| ) can be written 



Q. = 



H = 3H 2 



1 



R - A. 



(2.14) 



F or models with 3 R > 0, it is obvious that D = 
H 2 + 3 R/6 is a dominant quantity. We can thus use 
D, rather than H, to obtain compact variables for the 
k = 1 case: 



H 



A 
3~D 
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(2.15) 



Note that the sign of Q tells if a model is in an expanding 
or a contracting epoch. Introducing a new time variable 
defined by ' = D^ 1 d/dt, the equations become 



D' = -^-(l-n A )QD, 



Of 
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(i-n A ) 



(i-Q 2 ), 



f2 A = 3 7 (1 - «a)Q^a- 



(2.16) 

(2.17) 
(2.18) 



The evolution equation for D decouples, so that we ob- 
tain a reduced system in terms of Q and fl A . The equilib- 
rium points F and dS appear in this system as well. They 
correspond to (Q, Q\) — (±1,0) and (Q, Oa) = (±1, 1), 
respectively. Note that there is one expanding and one 
contracting version of each model. For 7 > 2/3 there is 
an additional equilibrium point E, given by 



Q = o, n A 



37-2 



(2.19) 



and with saddle stability. This is the Einstein static so- 
lution, for which S = 0, S = 0. This implies H = 0, 
H = 0, ft — ► 00 and q — > 00, see Pj. For each value S c of 
the scale factor, there is a value A c of the cosmological 
constant that results in an Einstein static universe. More 
precisely, 



Ar.= 



37-2 k 



(2.20) 



so in principle, the equilibrium point E represents a one- 
parameter set of Einstein static universes. 

The full dynamical system is obtained by matching the 
K > state space with one expanding and one contract- 
ing K < state space, see Figs. |l| and |^. The left half of 
the state space corresponds to expanding models, while 
the right half contain collapsing models. This is indicated 
by the subscripts of the various equilibrium points. Com- 
pact state spaces for all the Friedmann-Lemaitre models 
were first presented in , and also in |l^] . 
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FIG. 1. State space for Friedmann-Lemaitre models with 
2/3 < 7 < 2. In the K < regions (triangular), the vertical 
axis corresponds to Oa, and the horizontal axis to K. In the 
K > region (rectangular), the vertical axis corresponds to 
£7a, and the horizontal axis to Q. Subscripts on equilibrium 
points refer to the sign of H there. 

The state space looks different, depending on whether 
7 is less than or greater than 2/3. 

Referring to Fig. |l|, we comment on the different 
types of solutions for 7 > 2/3. The physical domain 
is bounded below by the invariant submanifold Q\ = 
and on the other sides by the invariant vacuum subman- 
ifold fl = 0. The bottom two apexes are the Milne uni- 
verse (Minkowski space-time in expanding coordinates) 
and the top ones are the de Sitter universe. For K < 
and H > 0, all solutions with O > and A > expand 
from an initial big bang singularity + F and evolve to the 
de Sitter model; this region is bounded by the K = in- 
variant submanifold (the straight line from + F to +dS). 
The time reverse region occurs for H < 0. 

For K > and H > 0, the situation is more compli- 
cated [|l9| . Models that start out with a sufficiently small 
value of A are closed Friedmann-Lemaitre models. They 
enter the contracting state space and recollapse to a 'big 
crunch' at _F. Models with large enough A will evolve 
to a de Sitter model. These are known as Lemaitre mod- 
els. In between these two classes, there is a separatrix 
orbit, corresponding to a model with A = A c and future 
asymptotic to the Einstein static universe. There are also 
models (starting off in the region with H < and crossing 
to H > 0) that contract and expand again, without con- 
taining a singularity. The separatrix between these and 
the Lemaitre models is the Eddington-Lemaitre model, 
which is past asymptotic to the Einstein static universe. 
The vacuum orbit from _dS to + dS corresponds to the 
de Sitter universe in different slicings. 

£1 = 
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FIG. 2. State space for Friedmann-Lemaitre models with 
< 7 < 2/3. In the K < regions (triangular), the vertical 
axis corresponds to 57a, and the horizontal axis to K. In the 
K > region (rectangular), the vertical axis corresponds to 
Oa, and the horizontal axis to Q. Subscripts on equilibrium 
points refer to the sign of H there. 

For 7 < 2/3, the situation is quite different. There 
is no longer an equilibrium point corresponding to the 
Einstein static universe. Also, the stability of the Fried- 
mann and Milne points have changed. For K < 0, all 
solutions are asymptotic to the Milne universe in their 
past, and evolve to the de Sitter model in the future. For 
K > 0, there are only singularity- free contracting and 
re-expanding models. 



B. The (q, SI) diagram 

A useful way to illustrate cosmological models is to 
construct a plot of the density parameter Q against the 
deceleration parameter q. This was first done by Stabell 
& Refsdal for the dust (7 = 1) Friedmann-Lemaitre cos- 
mologies Q, and extended to Friedmann-Lemaitre mod- 
els with matter and radiation by Ehlers & Rindlcr 
and to general equation of state by Madsen et al. 

In a (g, O) diagram, the different submanifolds of the 
Friedmann-Lemaitre models are represented by straight 
lines as follows: 



vacuum 
K = 

n A = 



n = 



q < 
-1<<?<^ 



37-2 



^1 



We present the diagram for Friedmann-Lemaitre models 
with 7 = 1. Diagrams for other values of 7 are given in 




2 The parameter that used to be called a |lj is the same as 

n/2. 



FIG. 3. (q, fi) diagram for Friedmann-Lemaitre models 
with 7=1. 

As the Friedmann-Lemaitre models constitute a two- 
dimensional system, the (q, fl) diagram will just be an 
alternative state space (but only representing part of the 
full space represented in Fig. |I]). This enables us to 
uniquely identify equilibrium points in this diagram, in 
correspondence to those in Fig. |l|. More e xplici tly, the 
evolution equation for Q is given by Eq. (2.13), while 
that for q can be expressed as 



q' = 2q(l 



.g)-^(3 7 -2)0. 



(2.21) 



For models represented by higher-dimensional dynamical 
systems, this one-to-one correspondence will in general 
no longer exist. The (g, f2) diagram then represents a 
two-dimensional projection of the full state space, and 
evolution curves may consequently cross each other. 



III. BIANCHI COSMOLOGIES 

By dropping the assumption of isotropy, more general 
classes of models are obtained. The Bianchi cosmologies 
are models that have a three-dimensional isometry group 
G3 acting on the spatial surfaces of homogeneity. There is 
a locally rotationally symmetric (LRS) subclass of these 
models, corresponding to a G4 isometry group. 

The anisotropy of these models manifests itself in a 
non- vanishing shear tensor a a p. In classifying the mod- 
els, one usually considers the objects n a p and a Q , which 
describe the structure constants of the symmetry group 
G3, see e.g. @. The evolution equations can be obtained 
e.g. , by appropriately specializing the equations given in 
[pifl . Upon performing an expansion normalization in 
analogy with Sec. ||, and with additional definitions 



— ^ 

= H ' 



N af 3 



- H ' 



A a = 



H 



(3.1) 



we obtain 



The Friedmann equation: 

f2 = l-£ 2 + K-fl A , 
The Raychaudhuri equation: 



(3.2) 



q = 2S 2 



37- 2 
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where 



S 2 = > 0, 



K = -A a A a - i [2N a0 N a ? - (N a a f] 



(3.3) 

(3.4) 
(3.5) 
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In addition to these equations, there are also evolution 
equations for £ Q/ g, etc. Here, K is the curvature parame- 
ter defined in Sec. |l[ The only Bianchi models for which 
K may be positive (i.e. 3 R > 0) are of type IX. For all 
other Bianchi models, it follows that f2, S 2 , —K, and fl\ 
are compact and only take values between and 1. In 
addition, — 1 < q < 2. Note that this does not imply that 
dynamical quantities like N a p need be compact. 

The fluid motion need no longer be orthogonal to the 
surfaces of homogeneity, i.e. Bianchi models may be 
tilted. In this case, the rest spaces of an observer co- 
moving with the fluid need not be homogeneous. When 
following the normal congruence, on the other hand, the 
fluid will no longer look perfect. Throughout this pa- 
per, kinematical quantities with respect to the normal 
congruence are used. A formalism for studying tilted ho- 
mogeneous cosmological models has been given by King 
fe Ellis @. 

In the orthogonal (non-tilted) case, the energy conser- 
vation equation takes the s ame form as in the Friedmann- 
Lemaitre case, Eq. ( [2.13 ), ensuring that those models 
have an invariant vacuum (f2 = 0) submanifold. In the 
tilted case too this will be an invariant sub-manifold (for 
normal fluids), as can be seen on using a comoving de- 
scription. 

Wald H has shown that Bianchi models with K < 
will asymptotically approach a de Sitter universe in the 
sense that they isotropize (S 2 — > 0) and A dominates. 
However, as pointed out by Raychaudhuri & Modak p3[ , 
it is not guaranteed that a tilted fluid will become par- 
allel with t he no rmal congruence in this limit. Indeed, 
in Subsec. Ill D we will present a counter-example for 



which the tilt is non-vanishing for certain equations of 
state. 



A. The (g, fi) diagram 

As discussed above, SI and q are compact for non-type- 
IX models. This is useful as it leads to a compact (q, f2) 
diagram, even though the state space (which usually in- 
volves more variables) may be non-compact. The bound- 
ary of the physical part of the {q, J!) diagram can be in- 
vestigated by pair- wise setting of S 2 , K and £l A to zero: 



£ 2 = =K 
£ 2 = =Q A 

k= o =n A 



_2 

"=3^2? 
^ = 3(2-7) ( 2 



The resulting diagram is a triangle with its apex at 
for 7 = 



- ^7-2 



1. 



1). In Fig. |, this region is depicted 



An important consequence of the compactness of the 
(q, ft) diagram for non-type-IX models is that equilib- 
rium points with non-zero fl A (which, by Eq. ( pT2| ), 
must have q = —1) necessarily have (q, ft) = (—1, 0). In 
addition, it follows that this uniquely specifies the vari- 
ables (£ 2 ,AT,O0 to be (0,0,1). 




-l o i 2 q 

FIG. 4. Physical domain of the (q, Q) diagram for Bianchi 
models with K < and 7 = 1. 

In the case of orthogonal Bianchi models (including 
type IX), it is possible to show that the de Sitter point 
is the only new equilibrium point when extending these 
models from the A = case (with the exception of possi- 
ble equilibrium points at infinity for type IX), and that 
this point always is a sink. This follows from analyzing 
the general evolution equations of the orthogonal models, 
given in e.g. Jig ]. As we will see, tilted models may have 
additional equilibrium points. 



B. Orthogonal type I 

The simplest Bianchi models are the orthogonal mod- 
els of type I. The shear tensor is diagonal, and conse- 
quently there are only two independent components. In 
expansion-normalized form, we take them to be 



£+ — 2^ 22 



£33), £- = 



1 



2 v / 3 



(£ 



22 



J 33j 



(3.6) 



With these definitions, it follows that S 2 = £ 



Furthermore, N a p = = A a , so that K = 
duced dynamical system becomes (see e.g. 

(B) 



SI a 



-(2-g)E±, 
2(l + q)Q A . 



0. The re- 
^ and Eq. 

(3.7) 
(3.8) 



This three-dimensional dynamical system is compact, 
and has the following invariant submanifolds: 

S+ = (722 = — O33 

E_ = the LRS submanifold 
Q = the vacuum boundary 
SIa = the A = submanifold 

Of these, the last two constitute the boundary of the state 
space. Apart from the equilibrium points in the A = 
submanifold (see, e.g. |l(|), the de Sitter (dS) point is 
the only additional equilibrium point, as expected. 
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F 
K 
dS 



Fricdmann 
Kasner circle 
de Sitter 



saddle 
source 
sink 



The state space is depicted in Fig. ||. 

dS 




K 

FIG. 5. State space for orthogonal Bianchi models of type 
I. The vertical axis corresponds to Qa, and the horizontal 
planes are spanned by £+ and E_. The 'bottom' is the 
Qa = submanifold. 

As the dynamical system is invariant under rotations 
around the F-dS axis, there will be a one-to-one cor- 
respondence between the state space and the (q, Q) 
diagram (Fig. ||), even though the system is three- 
dimensional. Note that the Kasner circle is represented 
by a single point in that diagram. 




-l o i 2 q 

FIG. 6. The (q, Q) diagram for orthogonal Bianchi models 
of type I with 7 = 1. 



C. Orthogonal type II 

This model also has the shear fully specified by £+ 
and £_, defined in Eq. (|]^). In addition A a = 0, and 
there is only one non-zero component of N a g. W e will 
take it to be N = N u > 0. It follows from Eq. that 
K = —N 2 /12. The redu ced d ynamical system becomes 
(see e.g. |J|tO| and Eq. (|t|)) 



(3.10) 
(3.11) 
(3.12) 

The in- 



£-' = -(2 -<?)£_, 
N' = (q— 4£+) A, 
Q A ' = 2(l + q)n A . 

Again, we end up with a compact state space, 
variant submanifolds are the following: 

N = the type I submanifold 

£_ = the LRS submanifold 

O = the vacuum boundary 

51a = the A = submanifold 



There is one new equilibrium point when 7 > 2/3: P±~ , 
described in e.g. pi]. Although a sink in the A = 
submanifold, the P{ point is a saddle in the full state 
space. The future attractor of the system still is the dS 
point, in correspondence with the discussion of this point 
in Subsec. 



[II A 



As the state space is four-dimensional 
and hard to visualize, we present a state space diagram 
for the LRS submanifold (Fig. fj]). Note that this is a 
special case of the two-fluid models studied in 0] . 




K F -+ 

FIG. 7. State space for orthogonal LRS Bianchi models of 
type II. The vertical axis corresponds to Qa, the line between 
K_ and K+ is the £+ axis, and the third direction is the iV 
axis. The 'bottom' is the Qa = submanifold. 




£+' = -(2-g)£+ + iv 2 A 



(3.9) 



-1 1 2 q 

FIG. 8. The (q, Q) diagram for orthogonal LRS Bianchi 
models of type II with 7 = 1. 

To obtain the corresponding (q, Q) diagram, we have 
studied the orthogonal LRS type II system numerically. 
Some sample curves are displayed in Fig. ||. Note that 
the equilibrium points still can be uniquely identified, 
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except the two Kasner points, which both correspond to 
the point (2, 0). The curves that have a cusp or loop 
correspond to state space orbits that start at K_ and 
approach K + , from which they are subsequently repelled. 
The vacuum orbit between K_ and K + has q going from 
2 to and back to 2 again, hence the arrows in both 
directions on the positive part of the g-axis. 



D. Tilted LRS type V 

So far, we have only provided examples of models 
where the fluid flow is orthogonal to the surfaces of ho- 
mogeneity. However, as discussed above, Bianchi models 
allow for a fluid flow that is tilted with respect to the nor- 
mal congruence of the symmetry surfaces. Here, we will 
consider tilted LRS type V models. The A = case of 
these models were studied in J2j| . Our starting point will 
be the dynamical systems analysis of tilted type V mod- 
els by Hewitt & Wainwright G(|, which we will extend 
to the A^O case. 

As previously pointed out, the kincmatical quantities 
associated with the normal congruence n a of the spatial 
symmetry surfaces, rather than the fluid flow u a , are used 
as variables. Following King & Ellis ^2|, a tilt variable 
v is introduced, so that in an orthonormal frame where 
n a = (1,0,0,0), we have 



1 



(M,o,o). 



(3.13) 



In the LRS case, the normal congruence is fully described 
by its expansion 9 and one shear variable, taken to be 
a = 3[(T22 + C33]/2. In addition, the structure constants 
are fully described by one variable, which is taken to 
be a\. The field equations and conservation equations 
now give a constrained dynamical system in the vari- 
ables (9,(T,ai,v). The relations between (9, a) and the 
kinematic quantities of the fluid flow have been given in 
p6| . By introducing dimensionless variables E = a/9 
and A — 3ai/#, and a new dimensionless time variable 
' = 3/9d/dt, the 9 equation decouples, and a reduced 
dynamical system with compact state space is obtained 
|26| . A cosmological constant is included by introducing 
an additional variable fl\ = 3 A/9 2 . The equations be- 
come 



The Friedmann equation: 



3^» 
9 2 



= 1 - E 2 - A 2 - A , 



(3.14) 



where \i n is the energy density as seen by an observer 
moving with the normal congruence. 



The Raychaudhuri equation: 
e' = -(l + q)9, 



(3.15) 



9 9.2 o 3(2- 7 ) + (5 7 -6)« 2 
q = 2 - 2A 3n A - 2[1 + (7 _ 1>2] A, (3.16) 



Reduced system: 

E' = -(2- g-2Au)£ 
A' = (q + 2Z)A, 
v(l-v 2 ) 



(3.17) 
(3.18) 



l-(7-l)v 2 
fV = 2(1 + q)Q A , 
Constraint: 



2E + 37-4- 2(7- lL4u], (3.19) 

(3.20) 



jvtt + 2[l + (j-l)v 2 ]Ai: = 0, (3.21) 

There is also an auxiliary equation for Q' that shows 
that = is an invariant submanifold. Note that the 
curvature variable used elsewhere in this paper is given 
by K — —A 2 . We assume that A > 0, without loss of 
generality. 

The system has a number of invariant submanifolds: 

S = CT22 = — 033 

A = the LRS type I submanifold 

v — the orthogonal submanifold 

v = ±1 extreme tilt 

Q = the vacuum boundary 

£1a = the A = submanifold 

Note that v — implies that either E = or A = 0. 
This follows from the constraint, Eq. ( p. 21 ). The former 



case results in the same equations as for the Friedmann- 
Lemaitre models with K < 0. 

The following equilibrium points can be identified: 
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n A 
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flat Friedmann 
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Milne 





1 
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±1 





M 
K± 


!J 

Kasner 




±1 


1 




37-4 

2(7-1) 








KI 




±1 





Tl 





dS 


de Sitter 











1 


dS± 


1; 








±1 


1 



Note that all these equilibrium points except F are on 
the vacuum boundary (This does not mean that fj, need 
to be vanishing. Indeed, the Kasner points have /i — ► 00, 
9 — > 00). In addition, there is a line H, of equilibrium 
points at v = 1 for which A = 1 + E, -1 < E < 0. 
Furthermore, there are two equilibrium points K^ for 
which E = v = ±1. As they are not part of the closure 
of the interior of state space, we need not consider them. 

Doing a stability analysis, most of the equilibrium 
points are found to be saddles. The exceptions are the 
following: 

0< 7 <2/3 2/3< 7 <4/3 4/3 < 7 < 2 
M + source source source 

K+ saddle source source 

K7 source saddle saddle 

Kl source source source 

dS sink sink saddle 

dS ± saddle saddle sink 

H source source source 
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Thus, it is found that there are three equilibrium points 
for which £l\ ^ 0: the de Sitter point dS, and two de- 
Sitter points dS^ with extreme tilt (v — ±1)- The vac- 
uum orbits between these points and dS correspond to 
the de Sitter model in different slicings, i.e. with differ- 
ent values of v. When 7 > 4/3, it turns out that dS is 
a saddle point, while dS^ become sinks. The threshold 
value 7 = 4/3 is particularly interesting, as = 1 then 
implies v' — for any value of v. In other words, for each 
value of the tilt variable v, there is a de Sitter equilibrium 
point. Figs. |9|— [l~T| depict the state space of these models 
for various 7-intervals. As this state space represents a 
dynamical system with four variables and a constraint, it 
is not always possible to unambiguously identify a direc- 
tion with a certain variable. Still, the vertical direction 
corresponds to VIa- 



ds + dS" 




FIG. 10. State space for tilted LRS Bianchi models of type 
V with 7 = 4/3. The vertical axis corresponds to Oa- The 
'bottom' is the Qa = submanifold. 



dS + dS" 




K F K + 

FIG. 11. State space for tilted LRS Bianchi models of type 
V with 4/3 < 7 < 2. The vertical axis corresponds to Oa- 
The 'bottom' is the Qa = submanifold. 



Thus generically the space-time becomes de Sitter-like, 
in accordance with the Wald theorem, but the tilt does 
not die away, in accordance with the cautionary note 
given by Raychaudhuri and Modak; hence at first glance, 
isotropization of the cosmology does not occur. However 
one must be cautious about this: de Sitter space-time can 
be described in terms of many congruences of curves, be- 
cause there is no preferred congruence in a space-time of 
constant curvature (cf. ^7j); on switching to a descrip- 
tion based on the family of spacelike surfaces orthogonal 
to the flow lines (highly tilted relative to the original sur- 
faces of constant time), the space-time may again appear 
isotropic at late times. But if so, that is not the end of 
the story: the question then occurs as to what happens 
on reheating at the end of any inflationary epoch that 
may be represented by this model: does the matter in 
the universe remember the original reference frame? If 
it does, then after reheating the universe may be highly 
anisotropic; but if it does not, then the matter may still 
emerge from reheating in an anisotropic way precisely be- 
cause there is no preferred reference frame in a de Sitter 
space-time, so the matter does not know what reference 
frame to choose after reheating (sec |28j for a discussion) . 

A further complication is that it may be possible to 
extend beyond the apparent limiting point to a region 
with timelike symmetries, as happens in the solutions 
discussed by Collins and Ellis (note that this is not a 
problem for the 7 = 4/3 de Sitter- like solutions with 
non-extreme tilt); this needs careful investigation.^] 

Thus the implications are unclear; they certainly de- 
serve further investigation, in particular because tilt is 
the generic case. 



3 The flow lines cannot continue across a surface Q, = 0, 
as this is an invariant surface of the flow, but one may be 
able to continue the family of solutions beyond by analytic 
continuation. 
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IV. KANTOWSKI-SACHS COSMOLOGIES 

In order to make the discussion of spatially homoge- 
neous cosmologies complete, we examine the Kantowski- 
Sachs cosmologies, which are models with an isometry 
group G4 whose G3 subgroup acts multiply transitively 
on two-dimensional spherically symmetric surfaces. The 
Einstein tensor is diagonal, showing that these models 
cannot be tilted. The global structure of these models 
has been studied by Collins [p9| , and a dynamical sys- 
tem with compact state space has been obtained |50],[n| . 
Here, we will extend this dynamical system to include a 
cosmological constant. 

The Friedmann equation for Kantowski-Sachs models 
can be written 



fi = 3H 2 - 3ct+ 2 + 3 R/2 - A, 



(4.1) 



where 3 R > and a + = (er 2 2 + 033 )/2. Following 
we assume /1 > and identify the dominant quantity 
\JH 2 + 3 R/6. We are now able to obtain a compact 



D 



state space by normalizing with D, rather than H . Thus, 
we define compact variables Qo = H/D, Q+ = <j + / D and 
n A = A/ 3D 2 . In addition, a new dimensionless time ' = 
D^ 1 d/dt is introduced. The dynamical system becomes 
The Friedmann equation: 



D 



= 1 - 



-n A , 



30 



n, 



3D 2 

Decoupled equation: 
D' = - Q+ + Q Q - Q Q+ 

Reduced system: 

Oo' = (l-Qo 2 ) (i + Q Q Q+-ZQ^ 

Q+ = -(1 - Q+ 2 ) (1 - Qo 2 + 3Q Q+) 
37 

+—^dQoQ+, 
D' - 

o A = -2— n A , 

In addition, we also obtain fin' = — (37Q0 + 2%-)Qd, 
showing Qjj = to be an invariant submanifold. Some 
other invariant submanifolds are also found. To summa- 
rize, we have: 



(4.2) 

D, (4.3) 

(4.4) 

(4.5) 
(4.6) 



2o=±i 
n D = 

Ov = 



flat ( 3 R = 0) submanifolds 
the vacuum boundary 
the A = submanifold 



The system has a number of equilibrium points: 

Qo Q+ stability 



±F 


flat Friedmann 


±1 








saddle 


+K± 


Kasner 


1 


±1 





source 


_K± 


)5 


-1 


±1 





sink 


+dS 


de Sitter 


1 





1 


sink 


_dS 


55 


-1 





1 


source 


±x 


see text 


±1/2 Tl/2 3/4 


saddle 



All these are points on the vacuum boundary, except ±F 
for which flu = 1. It should be pointed out that although 
±F are saddles in the whole interval < 7 < 2, the 
stability is changed at 7 = 2/3. We will only present 
state space diagrams for the case 2/3 < 7 < 2. 

The +X point corresponds to an exact solution with 
H = —cr+ = ^VA, and line element 



ds 2 



-dt 2 



e 2VXt dx 2 



dy + dz 2 



(4.7) 



This can be interpreted as a space-time where the shear 
exactly counter-balances the expansion in all but one spa- 
tial direction. This solution has a 'pancake' singularity 
in the past. The contracting analogue of this solution 
is associated with _X. The solutions corresponding to 
these points have been discussed in p2|-|35[ . 

Thus, the state space (Fig. |l2|) exhibits many similar- 
ities with the K > Friedmann-Lemaitre models: the 
state space is divided into two halves, one of which cor- 
responds to an expanding epoch (Q > 0) an d one where 
the models are contracting (Qo < 0). In each half, there 
is a point corresponding to the de Sitter solution. The 
points ±X have a role similar to the Einstein static point 
in the Friedmann-Lemaitre state space. Indeed, just as in 
that case, the Kantowski-Sachs state space also contains 
models that contract and expand again without ever be- 
ing singular. Also, note the surfaces of separatrix orbits 
associated with the Friedmann points ±F. 




K 

+ - . . 

FIG. 12. State space for Kantowski-Sachs models with 
2/3 < 7 < 2. The points ±X are located on the vacuum 
boundary. The vertical axis corresponds to JIa, the line from 
+ K_ to _K_ is the Qo axis, and the third axis corresponds 
to Q+. The 'bottom' is the (Ia = submanifold. The 
equilibrium points that have been drawn in shaded color are 
'screened'. 

To facilitate an easier understanding of the Kantowski- 
Sachs state space (Fig. |l2 | ), we present a picture of the 
f^A = submanifold (Fig. [l3[ ) . This state space was first 
given in 13^]. 







+ K + H>0 H<0 K + 




FIG. 13. State space for the Kantowski-Sachs Qa = sub- 
manifold with 2/3 < 7 < 2. The horizontal axis corresponds 
to Qo, and the vertical axis to Q+. 

V. CONCLUSIONS 

We have examined the evolution of various homoge- 
neous cosmological models. The inclusion of a cosmolog- 
ical constant greatly affects the late-time behavior. For 
example, most Bianchi models will evolve to a de Sit- 
ter universe and isotropize, as pointed out by Wald ||. 
The models presented here confirm the broad picture pre- 
sented by Wald of a positive cosmological constant lead- 
ing to isotropization in the orthogonal case, but leave 
it open in the more general tilted case. More precisely, 
there are models for which the tilt does not vanish at late 
times. To an observer moving with the fluid, this model 
will not seem to isotropize. These models may appear 
isotropic in another frame at late times, but whether that 
is in fact true or not is an open question. Tilted models 
deserve more investigation, as does the case of a realistic 
scalar field representation (the investigation here can rep- 
resent a model including a slow-rolling scalar field, but 
not one where the dynamics of the scalar field is more 
interesting) . 

ACKNOWLEDGMENTS 

The authors would like to thank Alan Coley, Ulf Nils- 
son, Claes Uggla and John Wainwright for helpful discus- 
sions. MG was supported by a grant from Wallenbergs- 
stiftelsens jubileumsfond. 



[1] R. Stabell and S. Refsdal, Mon. Not. R. Astron. Soc. 132, 
279 (1966). 

[2] M. S. Madsen and G. F. R. Ellis, Mon. Not. R. Astron. 
Soc. 234, 67 (1988). 



[3] M. S. Madsen, J. P. Mimoso, J. A. Butcher, and G. F. 

R. Ellis, Phys. Rev. D 46, 1399 (1992). 
[4] A. A. Coley and J. Wainwright, Class. Quant. Grav. 9, 

651 (1992). 

[5] R. M. Wald, Phys. Rev. D 28, 2118 (1983). 

[6] U. Brauer, A. Rendall, and O. Reula, Class. Quant. Grav. 

11, 2283 (1994). 
[7] C. W. Misner, Phys. Rev. Lett. 22, 1071 (1969). 
[8] J. Wainwright, G. F. R. Ellis, and M. Hancock, Class. 

Quant. Grav. 15, 331 (1998). 
[9] A. Rothman and G. F. R. Ellis, Phys. Lett. B 180, 19 

(1986). 

[10] J. Wainwright and G. F. R. Ellis, Dynamical systems 
in cosmology (Cambridge University Press, Cambridge, 
1997). 

[11] V. A. Belinskii, L. P. Grischuk, Ya. B. Zel'dovich, and I. 
M. Khalatnikov, Sov. Phys. JETP 62, 195 (1985). 

[12] V. A. Belinskii, H. Ishihara, I. M. Khalatnikov, and H. 
Sato, Prog. Theoret. Phys. 79, 676 (1988). 

[13] E. J. Copeland, A. R. Liddle, and D. Wands, Phys. Rev. 
D 57, 4686 (1998). 

[14] R. J. van den Hoogen, A. A. Coley, and J. Ibanez, Phys. 
Rev. D 55, 5215 (1997). 

[15] A. A. Coley, J. Ibanez, and R. J. van den Hoogen, J. 
Math. Phys. 38, 5256 (1997). 

[16] A. P. Billyard, A. A. Coley, and R. J. van den Hoogen, 
Phys. Rev. D 58, 123501 (1998). 

[17] S. W. Hawking and G. F. R. Ellis, The large scale struc- 
ture of space-time (Cambridge University Press, Cam- 
bridge, 1973). 

[18] J. Wainwright, "Relativistic Cosmology", in Proceed- 
ings of the 46th Scottish Universities Summer School in 
Physics, eds. G. Hall and J. Pulman, Institute of Physics 
Publishing, Bristol, 1996. 

[19] H. P. Robertson, Rev. Mod. Phys. 5, 62 (1933). 

[20] J. Ehlers and W. Rindler, Mon. Not. R. Astron. Soc. 238, 
503 (1989). 

[21] H. van Elst and C. Uggla, Class. Quant. Grav. 14, 2673 
(1997). 

[22] A. R. King and G. F. R. Ellis, Commun. Math. Phys. 
31, 209 (1973). 

[23] A. K. Raychaudhuri and B. Modak, Class. Quant. Grav. 

5, 225 (1988). 
[24] C. Uggla, Class. Quant. Grav. 6, 383 (1989). 
[25] C. B. Collins and G. F. R. Ellis, Phys. Rep. 56, 65 (1979). 
[26] C. G. Hewitt and J. Wainwright, Phys. Rev. D 46, 4242 

(1992). 

[27] E. Schrodinger, Space-Time Structure (Cambridge Uni- 
versity Press, Cambridge, 1961). 

[28] P. Aninos, R. A. Matzner, and T. Rothman, Phys. Rev. 
D 43, 3821 (1991). 

[29] C. B. Collins, J. Math. Phys. 18, 2116 (1977). 

[30] C. Uggla and H. von zur Miihlen, Class. Quant. Grav. 7, 
1365 (1990). 

[31] M. Goliath, U. S. Nilsson, and C. Uggla, Class. Quant. 

Grav. 15, 167 (1998). 
[32] E. Weber, J. Math. Phys. 25, 3279 (1984). 
[33] E. Weber, J. Math. Phys. 26, 1308 (1985). 
[34] 0. Gr0n and E. Eriksen, Phys. Lett. A 121, 217 (1987). 
[35] P. Vargas Moniz, Phys. Rev. D 47, 4315 (1993). 



10 



